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Abstract 

We consider the cosmological evolution of a brane in the presence of a bulk scalar field 
coupled to the Ricci scalar through a term f{(f))R. We derive the generalized Friedmann 
equation on the brane in the presence of arbitrary brane and bulk-matter, as well as 
the scalar field equation, allowing for a general scalar potential V{(j)). We focus on a 
quadratic form of the above non- minimal coupling —^((P'R and obtain a class of late-time 
solutions for the scale factor and the scalar field on the brane that exhibit accelerated 
expansion for a range of the non-minimal coupling parameter ^. 



1 Introduction 



In the last few years, an increasing number of cosmological data [1],[2| show that the 
universe undergoes accelerated expansion attributed to an energy component referred 
to as "dark energy". Dark energy is a significant percentage of the total energy of the 
universe. Although a cosmological constant is admittedly the simplest model of dark 
energy, the huge fine-tuning required for its magnitude makes theorists unhappy [3],[4j. 
Other attempted explanations of the origin of dark energy are phantom fields [S] and 
quintessence [6], as well as modifications of gravitational theory itself [7j. Independently 
of the dark energy puzzle, theories of extra spatial dimensions, in which the observed uni- 
verse is realized as a brane embedded in a higher dimensional spacetime (bulk), have re- 
ceived a lot of attention in the last few years. In this framework ordinary matter is trapped 
on the brane but gravitation propagates throughout the entire spacetime [8j.[9j.[lQj. 
The cosmological evolution on the brane is described by an effective Friedmann equa- 
tion [II],[l2] that incorporates non-trivially the effects of the bulk. Brane models provide 
us with new possibilities for the understanding of open cosmological issues like the ob- 
served accelerated expansion. The existence of a higher dimensional embedding space 
allows for the existence of bulk matter, which can certainly influence the cosmological 
evolution on the brane and be a major contributor to the dark energy. A particular form 
of bulk matter is a scalar field [13]. The presence of a bulk scalar field opens up the pos- 
sibility of a direct coupling of to the curvature scalar. A specific form of this coupling 
corresponds to the gravitational term appearing in the so-called tensor-scalar theory of 
gravity A bulk scalar field non- minimally coupled via a coupling of the form f{4>)R 
has also been considered in the Randall-Sundrum framework and a class of analytic as 
well as numerical solutions have been discussed [T5] . . [T7] . 

In the present article we study the cosmological evolution on a brane in the presence 
of a bulk scalar field, non-minimally coupled to the Ricci scalar 0. We derive the cosmo- 
logical evolution equation on the Brane as well as the corresponding scalar field evolution 
equation. In the latter we allow for a general scalar potential function. In addition to 
the scalar potential, these equations contain the unspecified function standing for 

the non-distributional part of the second derivative of the scalar field with respect to 
the fifth spatial coordinate. We focus on a quadratic choice of the coupling to the Ricci 
scalar, namely /(0) = 2M^(1 — ^0^/2). Assuming a simple quadratic form for the scalar 
potential on the brane and employing various ansatze for (j)"{t), we derive a class of ap- 
proximate late-time solutions. These solutions exhibit accelerated expansion for a range 
of the non-minimal coupling parameter ^. The plan of the paper is the following. In Sec- 
tion 2 we introduce the general framework of the model and derive Einstein's equations 
for arbitrary bulk and brane-matter. In Section 3 we calculate the evolution equations 
on the brane for a general coupling function. We focus on the specific case of a quadratic 
non- minimal coupling —^0^ R and write down the equations in a "dark energy variable" 
formulation. In Section 4 we discuss the late time behaviour in the absence of matter. 
Suitable choices for 0" and the scalar potential give accelerated expansion driven by a 
cosmological constant and a scalar field exponentially falling in time. Finally, in Sec- 

"'^A similar case with a 4-D scalar non-minimally coupled to induced Ricci curvature was treated in 
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tion 5, we present two approximate late time solutions for two different ansatze of (j)' 
exhibiting accelerated expansion for a limited range of the coupling parameter ^. 



2 General Framework 



Consider the Action 

S= j rf^a:v^|/(0)7^-A-l(V0f -V(0)-4")| 
+ / rf^Xv^|-a + 4-n, (1) 



describing gravitation in the presence of a scalar field that couples non-minimally 
to the bD Ricci scalar R through a general coupling function Standard Model 

matter, described by C'^^\ is confined on a Brane located at = y = 0. Additional 
bulk matter, distinct from the bulk scalar field (p, is described by C^^\ The 5D metric 
Gmn has signature (— , +, +, +, +). With f/^,^ we denote the iD metric on the Brane. 
Finally, a is the (positive) Brane-tension and A is the 5D cosmological constant. Varying 
the Action with respect to the metric we obtain Einstein's equations 

(7^Miv - IgmnTZ^ - VmV^/(0) + GMivVV(0) = ^Tmn , (2) 
where Tmn is the total energy- momentum tensor 

Tmn = T^MN + "^MN + "^MN " Gmn^ — G f^,y6'^j6'^a 6{y) . (3) 
'^MN stands for the scalar field part of Tmn, namely, 

Tj^l = ^m<PVn<P - Gmn ( ^ (V0)' + ) , (4) 

while the scalar field equation of motion is 

^9, dV ^df da ^, , „ 

- -7- + ^xr - 37'^(?/) = • (5) 



Let us now introduce the metric ansatz[ 

ds^ = -n^{y,t)dt^ + a'^{y,t)-fijdx'dx^ + b'^{y,t)dy^ (6) 

or 

f -nHy,t) \ 

Gmn = a\y,t)^,, ^ • (7) 

V b\y,t) J 

This metric corresponds to a Friedmann-Robertson- Walker geometry on the brane with 
a maximally symmetric 3-geometry 'jij. We assume Z2 symmetry and y takes values in 
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the [—00, +cxo] line. The functions n{y,t), a{y,t) and b{y,t) are continuous with respect 
to y but may have discontinuous first derivatives at the location of the brane. 

For the metric ansatz ([7]), the scalar field energy- momentum tensor is 



(<t>) 



/ 1^2 I I ^2 
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MN 
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where, we assume that the scalar field depends only on the fifth coordinate and denote 
by (j)' its derivative with respect to it and by its derivative with respect to time. The 
bulk and brane parts of the energy-momentum tensor can be parametrized as 



T 



MN 












(b) 



MN 















(9) 

where the hulk energy and momentum densities ps, Pb, Pb, as well as the energy ex- 
change function P5, are functions of time and of y, while the brane energy and momentum 
densities p, p are functions of time. 

Substituting the above ansatz and choosing Gauss normal coordinates {b{y,t) = 1), 
we obtain the equations of motion in the following form 



n~ \ h 

a 



a 



a 



+ ^0'' + Y^(0) + "^PB + y 5(1/)P + y + y A , (10) 
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^ 02 + 0'' + r(0) 
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2 2 

+ \l^APB - A) + \l^Ay) iP - ^(0)) > (13) 
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0+ (s---]^-n^S(P" + (!^ + 3-](P'\ +n^V' -n^nf' + n^a'S{y) = 0. (14) 
\ a n J y \n a J ) 

The Ricci scalar TZ appearing in the scalar equation is 

_ k 1 I a f ah I a" n" f a' n' 
7^ = 3— + — <^6- + 6 - -6 — } -6 2 6 - -6 



n^la \a j an \ a n V^/ ^ ^ 

We have denoted with dots the derivatives with respect to time and with primes the 
derivatives with respect to the fifth coordinate. Note though that V' denotes a derivative 
with respect to the scalar field 0. We have also allowed for a 0-dependence of the brane- 
tension a and denoted with a' its derivative with respect to the scalar field. 

Assuming Z2 symmetry and denoting a'{t) = a'{+0,t), n'{t) = n'{+0,t) and 
(f)'{t) = (j)'{+0,t), we may proceed to extract from the equations of motion the Junc- 
tion Conditions at ?/ = 0. Thus, from flTUl) . f[T5]) and f[T^ we obtain 

-6/--2f0' = i(p + a) 
a 2 

4/- + 2n7 + 2/V = ^(p-a) 
a 2 

-2(f)' + Af (n' + 3-^ = -a' 



Note that, from now on, /' stands for a derivative with respect to (j). We have also chosen 
?T.(0,t) = 1. These equations can be put in the form 

a' _ l (a + 2/V) l (3p-p) 1 (p + p) 

a 2 u 8 u 16 / ' ^ ^ 

/_ l (^ + 2/V) HSp-p) 3 jp + p) 

2 u 8 u 16 / ' ^ ^ 

u u 

where we have introduced 

u^2[^f + 8{ff] . (18) 
In the case that the brane tension does not depend on 0(t), these equations simplify to 

a' 1 (3p - p - 4(7) 1 (p + p) 
a 8 u lb J 

" 8 u + 16 / ' ^^"^ 

<^- = r^'^-^-'"^ (21) 
u 
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3 Evolution Equations on the Brane 



From equation ( ITTi) we may obtain on the brane a generalized continuity equation 



p + 3-(p + p) + 2P5 = 
a 



(22) 



that expresses energy conservation. Notice that this equation is /-independent. From 
equation (fT2l) we may obtain a second order or generalized Friedmann equation 



3/ 



+ 



+ 



k 



+ f + 3-/ + -(0)^--(y + A) = 

a 4 2 bAu 



6 {p + pf 1 /V 3/ 2 



(23) 



that expresses the time evolution of the scale factor a{t). We have taken k = 0. Finally, 
the scalar equation (fT^ gives on the brane 



+ 3 



+ 



- 0' U' + 3- 



a / \ a 



2/' 3- + n" 



(24) 



where 0"(t), 0'"{t), n"{t) stand for the non- distributional parts of these derivatives. The 
quantities a" and h" appear in equations ( fTOl) . (|T3l) considered on the brane and can be 
expressed in terms of 0" and standard quantities, i.e. a(t), 0(t), their time-derivatives 
and the various matter densities. Doing thafl we obtain the set of three independent 
equations, namely fl22l) and 



a a 
-+[- 



12/m 24/m dj 



1 d? 
+ _(y + A) _ J^(i + 4/") 



L 
3/ 



+ 3- 



(25) 



+ 3^0 - 0" + 1^02 ^ -^F'(0) - ^ (3P5 + Pb- Pb) 
a 2u u u 

+ 10— {V + A) + 4^^-^ ^ — ^— (p - 3p) + ^ ^ — ^— 



(26) 



We have enforced the low energy density approximation, where we neglect p^ terms. Note 
that the scalar evolution equation contains the unknown quantity 0". Although, up to 
now, we have worked in the framework of a general coupling function /(</>), we know 



^For simplicity we have taken <t' = 0. 
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that, since = this cannot be very far from a constant at late times (Einstein 

gravity). Thus, we may consider a quadratic couphng function 



= 2M3(^1-|02^ , (27) 

where ^ is a suitable coupling parameter |14] [T5]|16j . Note that fl27|) can be a good ap- 
proximation to a general coupling function /(0) ~ /(O) + /"(O)0^/2 + ■ ■ ■ for small 
<P « (2/(0)//"(0))i/2. 

In order to get the scale factor evolution equation (125|) into the familiar first order 
Friedmann equation form, we may introduce an auxiliary variable x{^)i ^ so-called dark 
energy field. In the dark energy field formulation, \l9l,'^^ equation (l25l) is replaced by a 
pair of two first order equations. These equations are 



, 2 

a 



27P + X + A, (28) 



X + A-<x + 



12/ 



Pb + 1(1 - 400' - - 2^0 U + 3-0 



= 47P5 - 27P - A . (29) 

Note that we have restricted ourselves to the (1271) form for the coupling function /, 
although we have not substituted this expression everywhere. The functions [3, 7 and A 
appearing in (|28l) and ( l29l) are defined as 



13 



1 



24/m 



, 7 = a/3, A 



12/ 



2u 



(30) 



The original second order differential equation fl25l) can be recovered by differentiating 
(1281) . inserting fl29l) and using the continuity equation fl22|) which is always assumed to 
hold. An equation expressing a/a in terms of the dark energy field can be written down, 
namely 



- = -lip + 3p + ^ 
a 6/ 



V - P 



1 



(1-400' + 2^0 + 3 



X + A. (31) 



The function A is the effective cosmological constant on the brane and, in the minimal 
case (/ = 2M^), the condition 



A = A + 



a 



24M3 







corresponds to the familiar Randall- Sundrum fine-tuning of no cosmological constant on 
the brane in the absence of matter. 
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4 Late Time Behaviour in the Absence of Matter 



Let us now consider the case of absence of any matter. In this case, our set of equations 
becomes 

'a' 



X + A, (32) 



a r 12/ 



1(1 _ 4^)02 _^ _2^0 U + 3^^ 



-A, (33) 



+ 3^0 - 0" + = -^y'(0) + 10^(F + A) + '^^^"+3^^">' . (34) 

Notice that in the minimal case (/' = 0) the scalar equation can be satisfied with a 
time-independent solution 0o of V^'(0o) = 4>" i provided 0" is time- independent as well. 
In this case, the dark energy equation (!33l) is just (Vq = V^(0o)) 

X + 4^fx-,JWl =0, (35) 



a V 24M3 
with a solution 



Finally, the effective Friedmann equation is 



aV ^ C Vo 1 cr^ 



We can always choose Vq = 0. Then, the condition for a vanishing cosmo logical constant 
on the brane takes the familiar Randall-Sundrum fine-tuning form A + (t^/24M^ = 0. 
In the case of vanishing cosmological constant the resulting dark radiation-dominaXed 
expansion is decelerating like a{t) oc t^/^. Actually a constant solution is possible in the 
presence of the non-minimal coupling as well and the above behaviour is not modified. 
Simply, the condition on the constant (f)" is replaced by the more complicated condition 
0" = QfV\(p)/u - 10/'(F + A)/n - 8/'(u + f'u')a^/u^. The dark energy field is now 
X = C/a^ + Vo/12/, while the Friedmann equation retains the same form with A defined 
in fl5Ul) . Depending on the fine-tuning condition that we impose, we either have a cosmo- 
logical constant dominated exponential expansion or a dark energy dominated a{t) oc t^^"^ 
expansion as in the minimal case. Since the effective cosmological constant A depends on 
both / and m, its value at late times will contain the coupling parameter ^. A residual 
cosmological constant will thus be present, even if we impose the Randall-Sundrum fine- 
tuning, due to the non-minimal coupling, as can be seen from the expanded Friedmann 
equation 



where we only kept terms up to 0^ and imposed the fine-tuning. Apparently, the presence 
of the Vq term is not necessary to get non-trivial results, except from the case of conformal 
coupling, i = ic = 3/32M3. 
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Instead of viewing the scalar equation flMl) as an equation that supphes us with a 
solution for for a given function 0", we may alternatively view it as an equation that 
determines the unknown function 0" in terms of a chosen configuration 4>{t). Thus, we 
may choose an exponentially decaying field configuration 



0(t) = 00 e 



-Kt 



(39) 



and assume that it is a solution of (IMI) for a suitable function 0". Substituting it into 
( 133!) . we obtain 



X + 4-<x + 



12/ 



— (1 - 400^ - V - 2^02 - 3k- 



-X. 



(40) 



or, using ( 132]) . 



A + 4- <^ X + 



24/ ' 



12/ 



+ 2^(X + A)0^ 



0. 



(41) 



At late times, we may expand the scalar potential in powers of ^ = (poe ""' as 

1/(0) ^ 1^(0) + \v"m' + . . . , 

assuming that it is a function of 0^. Thus, keeping only first order terms in 0^, we obtain 



X + 4- 



48M3 



-A(0) - 



+ /€2(l-8O-^"(0)-^n0) 



24M3 



,2— 



0, 



(42) 



where x = X + ^ and A(0) = (A + cr^/24M^)/24M^. To zeroth order, we have 

\/(0) 



x + 4- x-m 



24M3 







with a solution identical to (1361) . The resulting Friedmann equation is just (!37j) . Thus, 
at late times, depending on the fine-tuning condition imposed, we either have exponen- 
tial expansion driven by a cosmological constant or a dark radiation driven decelerated 
expansion. The non-minimal coupling does not play any role to this order. Corrections 
which depend on ^ occur once we include terms of order 0^. In any case, the existence 
of the exponential solution ( l39l) rests on the scalar equation and ultimately on 0". Ex- 
panding it in powers of and substituting, we get 



-0"(O) + 



- 3kH{0) - {0"(O)}i + V"iO) + ^ (ViO) + A) + 



+ 



-3nH'{0) - ^{0"(O)}2 



+ 



(43) 
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where we have introduced the Hubble parameter H = a/ a. We have introduced the 
expansions 

H ^ H{0) + H'{0)(f) + ... 

and 

0" ^ 0"(o) + {0"(o)}i0 + ^W\o)h<P' + .... 

Equation (03]) is satisfied for 0"(O) = 0, 

{4>"{o)h = ^'- + v"{o) + 1 {v{o) + A) + 

and 

{0"(O)}2 = -6«:i/'(0). 

If in addition we restrict the forms of 0" and V, we would be forced to fine-tune the 
parameters k, H{0) and H'{0). For 0" = and V = 0, equation fH5]) reduces to 



i/'(0)=0, if(0) = ^f«:2^^A + 
3k \ 3 



^ (44) 



9M3 

which provides an approximate solution of constant H, which depends on ^. 



5 A Class of Approximate Solutions 

Let us consider now our set of cosmological evolution equations in the presence of brane 
matter but in the absence of bulk matter, apart from the scalar field. For the late-time 
behaviour of the scalar field we shall adopt here a decreasing power law ansatz with an 
as yet unspecified power a, namely 

Ht) ^ ^ • (45) 

We may use the original second order scale factor differential equation instead of the 
dark energy field formulation since they are both equivalent. For the brane matter we 
shall assume the equation of state p = wp. With this equation of state, the continuity 
equation ( l22l) gives 

p + 3-(l + w)p = ^ p oc a-^(^+"') . (46) 
a 

For the above ansatz our pair of evolution equations gives 
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2u 



u 



+ 10— (y + A) + 4^^-^ ^ — '—p{l - 3u;) + ( — '— 



(48) 



For the chosen quadratic form of the couphng function fl7r|) . u takes on the simple form 
u = 12M3 [1 - ^ (1 - ^/Q 0V2J , where = S/32M^ is the conformal value of the 
non-minimal coupling parameter Q The coupling function derivatives that appear in the 
above equations are /' = — (2M'^)^0 and u' = — 12M^(1 — ^/^c)(,<P- We shall also assume 
a simple quadratic scalar potential and without loss of generality we shall take it to be 
vanishing at the origin, i.e. V{(f)) = /i^0^/2. 

Next, let us introduce an expanding scale factor ansatz 



a{t) 



(49) 



valid at late times, in terms of a power u to be determined. Inserting the scale factor 
ansatz fH9l) . the scalar field ansatz fH5|) and substituting the energy density P6l) as 

(50) 



into the Friedmann equation, we obtain 



v{2u 



48M3 6M3 



+ 4A + 



(t{1 -3w)CA 1 



12M3 



^3i/(l+w) 



+ 



cl 

f2a 



4^^4(1 -3w7) (2^0-0^ 1 



9 



2a 



-— [a{l - 16M^^) + 8M^^{3u - 1) ] 



(51) 



We assume that a > 1. Keeping only terms up to second order, we obtain the conditions 

2 



I' 



3(1 + m;) 
^2 



A + 







iy(2iy 



24M3 

a(l -3w;) 



C4 



(52) 

(53) 
(54) 



' 144(2M3)2 

Note that (|53|) is the standard fine-tuning for a vanishing cosmological constant on the 



brane. 

In a similar way the scalar equation gives 



a[a + 1 — Su] 



(3a+2)\ 



1 

2^ct^'' 



■^The effective cosmological "constant" on the brane is A = (A + o'^/2u) /12/. For the conformal 
value of the non- minimal coupling parameter £, — £,c, we have u = 12M'^ valid for all (p and the usual 
Randall-Sundrum fine-tuning gives A = as in the minimal case. 
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18M3 



t2 V 



2a 



We have used equation fl52l) . At this point we have to assume an ansatz for the unknown 
function 0". 

Ansatz 1 : (p" oc (p 

This amounts to a late time behaviour 



Substituting this into ( 155|) . we obtain the two relations 
(l-3«;)ea^ 



(56) 



Co 



-C1C4 = CiO; ( a + 1 — 3z/) , /i 



18M3 ' " ^ V ■ ^ ' A- 24M3 ■ 

We have considered terms only up to second order and ignored terms of 0(t~^") 
From equations flS^ . flMl) . flFTI) we may fix A, /i^ and C4 and obtain 



(57) 



3u - 1 



a 



± 



(3z/ 



(5^ 



Note that u is given by fl52|) . For z/ < 1, or equivalently, if > —1/3, the exponent a is real 
for any value of ^. In contrast, for u > 1, or equivalently, — 1 < < —1/3, the exponent 
a is real and positive for 



I2v{v - 1) 

This inequality restricts the values that ^ is allowed to take in order to have real values 
for the exponent a and a value of v corresponding to accelerated expansion. The presence 



of the free parameter jr- allows a certain degree of flexibility in the choice of ^ for small 



values of v. 

Ansatz 2 : 0" oc 

This amounts to a late time behaviour 



C2 



(59) 



Substituting this into the scalar equation (l55!l . we obtain 



1 



3 
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2a + ^(1-3^; 



2 

]_ 

t2" 



Neglecting terms of 0{t '^"), we obtain the relations 



9i 



2 I ^'^ A I 

+ —A + 



afa + 1 — Sv) 



3 ■ 9M3 ' 
(1 - 3w)^a 



(60) 



(61) 



(62) 



18M3 

From these relations we obtain an almost identical expression for the exponent a as in 
the case of the previous ansatz. The only difference is that the term C2/C1 is missing, 
namely 

(63) 



a 



± 



4 



3^z/(2z/-l) 



Note that v is given by ( l52l) . For z/ < 1, or equivalently, w > —1/3, the exponent a is real 
for any value of ^. In contrast, for i/ > 1, or equivalently, —1 < w < —1/3, the exponent 
a is real and positive for 

«<t|^. (64) 

llv[v — 1) 

In comparison to the previous ansatz, the range of allowed values for ^ is smaller here. 
In this case the free parameter ^ is present in the expression for jj? obtained from (!6T|) . 
namely 

^ Ci 24M3 ■ 



6 Conclusions 



In the present article we studied the cosmological evolution on a brane embedded in 
a bD bulk in the presence of a bulk scalar field non-minimally coupled to the Ricci 
curvature scalar through a term f{(f))R. We derived the cosmological evolution equations 
on the brane in the presence of arbitrary brane and bulk matter. The scalar field evolution 
equation on the brane contains the non-distributional part of the second derivative of the 
scalar field {4>")- The scalar evolution equation includes an unknown non-distributional 
part of the second derivative of the scalar field with respect to the fifth coordinate, a 
quantity that requires the knowledge of the dependence of the scalar field on the fifth 
coordinate. We have proceeded, considering this quantity as a function depending on 
the structure of the bulk and introduced various ansatze for it. Although we derived the 
evolution equations for a general coupling function, we focused on a quadratic form for 
it, namely /(</>) = 2M^{1 — ^(j)'^/2). We considered the "dark energy field" formulation of 
the scale factor evolution equation in which the second order differential equation for the 
scale factor is replaced by an effective first order Friedmann equation and an equation for 
an auxiliary variable. We first consider the case that no brane or bulk matter were present. 
We showed that late-time solutions exist, with a constant or exponentially decaying scalar 
field, in which the scale factor evolution is driven either by an effective cosmological 
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constant or by dark radiation. The value of the effective cosmological constant may in 
both cases depend on the couphng parameter. Although, no restriction is placed by 
these late-time solutions on the non-minimal coupling parameter, they do depend on the 
behaviour of 0". Next, we have considered the case of brane matter in addition to the 
scalar field. We introduced an expanding power law ansatz for the scale factor (a ~ t'^)and 
assumed a decreasing power law type of time-dependence of the scalar field at late times 
(0 ~ The resulting constraints depend on the assumptions on the behaviour of 0". 
Considering first the case 0" oc 0, we see that we indeed have such a late-time solution 
only for brane matter with an equation of state parameter —1 < w < —1/3 and a non- 
minimal coupling parameter ^ < .^c[(3z^ — 1)^ + 4C2/Ci]/12z/(z/ — 1). The ratio C2/C1 is 
the 0"/0 coefficient ratio. We next considered the case 0" oc and found that, again, 
such a late-time solution only for brane matter with an equation of state parameter 
—l<w<—l/3 and a non-minimal coupling parameter ^ < ^c(3z^ — l)^/12i^(i^ — 1). 
These results correspond to a simple quadratic scalar potential choice but are expected 
not to be modified by higher potential terms at late times. Nevertheless, they could be 
modified by the presence of bulk matter interacting with the brane. Depending on the 
ansatz for the bulk pressure and exchange densities, the range for w could be replaced 
by standard values. 
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